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NEW APPLICATIONS OF MACLAURIN’S SERIES 
IN THE SOLUTION OF EQUATIONS AND 
IN THE EXPANSION OF FUNCTIONS. 

BY P. A. LAMBERT. 

{Read April 3 , 1903.) 

I.—Introduction. 

The modern theory of differential equations is based on 
the expansion by Maclaurin’s series of the solutions of the 
equations in infinite series. The striking analogy existing 
between the theory of algebraic equations and the theory of 
differential equations suggested the possibility of expressing 
the solutions of algebraic equations in series to be obtained 
by an application of Maclaurin’s series. After some experi¬ 
menting the author happened on the device of introducing 
a factor.r into all the terms but two of the equation f (y) = 0, 
whereby y becomes an implicit function of x. The succes¬ 
sive ^-derivatives of y are now formed, and together with y 
are evaluated for -r = 0. By Maclaurin’s series the expan¬ 
sions of y in powers of x become known. If x be made 
unity in these expansions, the roots of f(y)~ 0 are found, 
provided the resulting series are convergent. 

To illustrate this method, consider the equation 

( 1 ) V * — 3?/ 2 + 75?/—10000 == 0 . 

Maclaurin’s series 


_ _i dy ° jl . i i — I 

V 2A) T ( i X(f x *r d,V(^ 2 ! ' dx 0 z B! ' dx 0 4 4 ! ' * * ‘ * 

S} W- Si •••• the ’lo¬ 
ot y, j i 0.when x is made zero, expands y, 

a function of x , in powers of x. 

By introducing a factor x in the second and third terms of 
(1) an equation is formed 

(2) y 4 — 3.r y 1 + 7 oxy — 10000 = 0 
which defines y as an implicit function of x. 


88 


LAMBERT—M ACL A l! It IN'S SERIES OF EQUATIONS. [April 3. 


Differentiating’ (2) twice in succession 
(3) 4 y 


dX 


8tf , +76y-teyg + 75*g = 0 




„ d-p . d~y A 
Gx^-^+To.r - = °. 


Making x zero in (2), (3) and (4) 

!/„=+ 10.— 10, + I0i /= T. —101 7 ^! 

^ = -.1125,—.2625, + . 1875 — .075 /ZZJ, + .1875 + .075/=! 

ax 0 

— —.0029, —.0029, —.0000015 + .0039/—1,— .0000015 —,0039i / :^I. 

«J 0 * 

Substituting these four sets of values in Maclaurin’s series 
and placing x = 1, the roots of equation (1) are found to be 

Pi 4~ 9.88b, P 2 == — 10.261, y$ = -f- .18*5 -j- 9.927 \f —1, 

Pi — - f- .1875 — 9.927 V—l, 

all correct to the last decimal. 

This method will be applied to the solution (II) of tri¬ 
nomial algebraic equations, (III) of general algebraic equa¬ 
tions, (IV) of trinomial transcendental equations, and finally 
(V) the method will be applied to obtain expansions com¬ 
monly obtained by Lagrange's series. 


II.— Trinomial Algebraic Equations. 

The general trinomial equation of degree n has the form 
(1) y n — nay u -* — 5 = 0. 

Introducing a factor x in the second term of (1) 

(2) y n — naxy u ~ k — 5 = 0. 

Applying the method and denoting the 11 th root of b by <0 
(3) y = o) + w l ~ k a -f- w I-2k (1 — 2k -f- n) 

+ (l _ 3 * + n) (1 — 3 k + 2;i) 

+ M i-**(1— 4/fc+n)(l — 4t + 2tt)(l— 4* + 3re)-^- + . . . . 
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To determine when series (3) is convergent, group the 
terms numbered 1, n + 1, 2?i + l, 3n + l, . . . . , then those 
numbered 2, ?7 + 2, 2?i + 2, 3?i + 2, . . . . , finally those num¬ 
bered n, 2 n, 371, 4?i,. Each of these n partial series 

is found by Cauchy's ratio test to be convergent when a n is 
numerically less than & _k ( 7 i—fc) k-n 6 k . When this condition 
of convergency is satisfied series (3), by substituting for u> 
in succession each of' the n values of the n lh root of b , deter¬ 
mines the n roots of equation (1). 

By introducing the factor x in the third term of equation 
(1) and applying the method a series is obtained which deter¬ 
mines k roots of equation (1), and by introducing the factor 
x in the first term of equation (1) a series is obtained which 
determines n — k roots of equation (1). The two series thus 
obtained are convergent when a n is numerically greater than 
k~*(n — &) k_n 6 k . When a n = k~*(n — &) k-n 6 k equation (1) has 
equal roots. There is therefore developed a complete theory 
of trinomial equations. 

The general fifth degree equation can, by Tschirnhausen 
transformations requiring the solution of equations of the 
second and third degrees only, be transformed into the tri¬ 
nomial equation y 5 + ay + b = 0. If a 5 is numerically less 

625 

than o^-6 4 , the five roots of this equation are found by apply¬ 
ing the method to \f + axy + 6 = 0. If a 5 is numerically greater 
than ||| 6 4 , the five roots are found by applying the method 

to y i +ay-\-bx = 0 and xy*-\-ay + b = 0. If a 5 numerically equals 
625 

0 ^6 4 , the fifth degree equation has equal roots, and the re¬ 
moval of the equal roots makes the solution of the fifth 
degree equation depend on the solution of an equation of a 
degree not higher than the third. A third degree equation 
becomes trinomial by removing the second term, which is 
accomplished by a linear transformation. The method of 
this paper therefore effects "the complete solution of the 
general fifth degree equation in infinite series. 

In Weber's Algebra, volume I, pages 396-399, the real and 
imaginary roots of the equation x 3 — 2x — 2 = 0 are computed 
by a method invented by Gauss for the solution of trinomial 



88 LAMBERT—MACLAURIN’S SERIES OF EQUATIONS. [Aprils, 

equations. The convergency test shows that the series found 
by introducing the variable factor in the second term is con¬ 
vergent. Now the mathematician is satisfied when the con- 
vergency of the infinite series he uses is established, but the 
computer desires that the infinite series he is obliged to use 
shall converge rapidly. By transforming the equation 
x 3 — 2x — 2 = 0 into another lacking the first power, which is 

3 

accomplished by placing x= ^the equation 

54y 3 — lSy — 23 = 0 is found. The series found by applying 
the method to 54y 3 —1 Sxy — 23 = 0 converges much more 
rapidly than the series obtained from the original equation. 
Differentiating 54 y 3 —IS xy — 23 = 0 four times in succession 
and making x zero, 

y 0 = .7524, — .8762 ± .3762/—3 

d p>-~ .1477,— .0738 =F.07381/_3 

dx o 


= _ .0019, + -0010 =f .OOIO/H 3 
— .0004, — .0002 T .00021/_3. 

4 ! dxj 

The three values of y are .SS9_and—.4492 ± .3012/—3, 
the corresponding values of x are 1.76S and 
— .8S47 “F .5898^ —1. If the computations are made by 
logarithms they are not very lengthy. 

The equation y 4 — 11727 7/4-40385 = 0 occurs in a paper by 
Mr. G. H. Darwin 4 ‘On the Precession of a Viscous Sphe¬ 
roid/’ published in the Philosophical Transactions of the Royal 
Society , Part II, 1879, page 508. The convergency test shows 
that the factor x must be introduced in the last and in the 
first terms. The equation therefore has two real positive 
and two imaginary roots. Applying the method to 
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y' — 11727)/ + 40385a - 0, 
y<s — 23.720,—11.360 ±11.360/^3 


T°«=— 1.148, — 1.148 


1 

2! 

<Py» 

.116, 

.058 dr 

1 

d 3 y 0 

.019, 

.010 =F 

3! 

(W 

1 

SL 

<*sT 

© 

.004, — 

.004 

4! 

1 

** 


.058/—3 
.010/^3 


Three roots of the equation are 21.432 and 
12.444 ± 19.759* 7 —1. Applying the method to 


xy i _11727)/ + 40385 = 0, 

i/ o = 3.4430, ^ = .0120,^ -^- = .0002. 
dx o (lx q 


The fourth root of the equation is 3.4558. 

This method applied to trinomial equations proves that an 
equation of degree n has n roots, determines how many roots 
are real, and presents a uniform scheme for computing all 
the roots, real and imaginary. 


III.— General Algebraic Equations. 

The method applied to the complete equation of degree n 
furnishes n series, and it becomes necessary to deter¬ 

mine which of these series give n convergent series for the 
roots of the equation and if possible to insure rapidity of 
convergence of these n series. 

Suppose the equation of degree n to be 

ay n + ftp y n ~ l + . . . . + flk y n ~*+ l + by D ~ k + b Q y n ~ k ~ l + .... 

+ 5 1 y n - k -^+ 1 + cy n ~ k ~ 1 + c 0 y n - k -'- 1 + . . . . + c m ?/ n - k - 1 - m + l 

+ dy n -*- ] ~ m + d Q y n - k ~ ] - m ~ l + . . . . + ry + s = 0 , 

and suppose the terms which are underscored to be the terms 
from which the two terms into which the factor x is not intro¬ 
duced must be selected by taking consecutive terms in regu¬ 
lar order from the left. The problem is how to recognize 
the terms which must be underscored. 
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If the factor x is omitted from the first two underscored terms 


i 

= ; if from the second and third underscored terms 


i 


*-K) : = 


if from the third and fourth underscored terms 
if from the last two underscored terms 

i 

_1_m Altogether n values of y 0 are found, and it 


is seen at a glance what values of y 0 are real and what are 
imaginary. In order that these values of y 0 shall be close 
approximations of the roots of the given equation, the suc¬ 
cessive derivatives d f-> d i ~, dj ~, • • • • must be small. 

dx 0 dx Q l iUq dx o 4 

i 

Forming Qh. corresponding to ?/ 0 = (— k and assuming that c 
dx 0 \ a* 


is of such a magnitude that the term containing c overshadows 

all the other terms in the numerator of it is found that 

dx o’ 

^ - is necessarily small if the ratio of b k+x to a'c k is numeri- 
dx n 


cally large. This same condition insures that the following 
derivatives 4^-, .... are small. 

UXq (IXq 

In like manner it is shown that the derivatives correspond- 
i 

( c \ 1 

ing to b ) are small provided the ratio of c x + m to 

b m d [ is numerically large, and that the derivatives corre- 


( d\ m 

—- ) are small provided the ratio of d n ~ ]k—l1 

to c m s n k i m numerically large. This ratio should, if 
possible, be made larger than 10 to insure rapid convergence. 

The directions for underscoring terms are therefore as 
follows : 

Underscore the first and last terms of the equation. Such 
other terms are to be underscored as satisfy the condition 
that if any three consecutive underscored terms be chosen, 
the ratio of the coefficient of the middle term with an ex- 
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ponent equal to the difference of the degrees of the first and 
third terms to the product of the coefficient of the first of 
the three terms with an exponent equal to the difference of 
the degrees of the second and third terms and the coefficient 
of the third term with an exponent equal to the difference 
of the degrees of the first and second terms shall be a large 
number. 

To illustrate the method, the following equations are 
discussed : 

(a) y b — 10 y 3 ■+- Gy + 1=0. 

Here all the terms are underscored, for the ratio of 
10 4 to 6 2 is large, and the ratio of 6 3 to 10 is large. The 
method must be applied to (1) y* — 10y 2 + 6xy + x = 0, 
{2) xy 5 — 10z/ 3 + 6?/ + x = 0 and (3) xif* — 10x*/ 3 + 6j/ + 1 = 0. 
The computation determines the following values: 


From 

a) 

Vo = 

+ 

3.167, 

— 3.167 



= 

dx Q 

— 

0.100, 

+ 0.090 



ii 

— 

0.008, 

+ 0.008 ; 

From 

(2) 

Vo = 

+ 

0.775, 

— 0.775 



dy 0 _ 
dx 0 

+ 

0.107, 

+ 0.060 



IO|H 

11 

— 

0.006, 

+ 0.016 ; 

From 

(3) 

Vo = 

+ 

0.166, 

i 

ii 


The roots of the given equation are ^ = + 3.05, 
y 2 =— 3.06, y 8 “ +0.S7, y 4 = —0.69, y 5 = - 0.17. 

( b ) x 4 -+ — 4.r 2 — llx -+• 4 = 0. 

Here the terms to be underscored in addition to the first and 
last are probably the second and fourth, but as the ratio of 
4 3 to 11 is rather small, it is safer to transform the equation 
into another lacking the second term by the substitution 
x = y — 1. There results 

y* _ 1Q£ + + 8 = 0. 

The terms to be underscored are the first, second and last 
and the roots are obtained by applying the method to 
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y x — 10?/ 2 + o.n/ + S.r = 0 and ,v 1 / — 10?/ 2 + oxy + 8 = 0. From 
each of the two equations two real roots, one positive and 
one negative, are found. 

(c) 7.c 4 + 20.r 3 + 3.c 2 — 16.c — 8 = 0. 

Here the terms to be underscored are probably the first, 
second and last, indicating the existence of two imaginary 
and two real roots, one positive and one negative. All doubt 
is removed by transforming by x = y — .7 into 

7y» + Aif — IS 42?/ 2 — 1.404?/ — . 509 3 = 0. 

The transformation x = y — .7 is selected because it is a 
simple transformation which makes the coefficient of the 
second term very small. 

(cl) x 5 + 12x 4 + 5 Or* + 150a; 2 + 201a; — 207 = 0. 

Here probably only the first and last terms are to be 
underscored, indicating the existence of four imaginary roots 
and one real positive root. Transforming by x = y — 2, 
which makes the coefficient of the second term small, 

2Z 5 . + W + 3/ + 4y 2 + by — 321 = 0. 

The roots are found by applying the method to 

y b + 2xy* + 3a;^ 3 + 4xy 2 + 5 xy — 321 = 0 
(e) x* — 80a; 3 + 199S.r 2 — 14937a; + 5000 = 0. 

Here probably every term should be underscored, indicat¬ 
ing four positive real roots. Transforming by the substitu¬ 
tion x =?/ + 20, 

y* _ 402 y 2 + 983 y + 25460 = 0. 

Here the terms to be underscored are the first, second and 
last. More rapidly convergent series are found by reversing 
the last equation, 

25460 c 4 + 9S3c‘- — 402 0 + 1 = 0, 

and making the substitution v = z — .01, whence 

2 5460s 4 — 35.4s 3 — 416. 214s 2 + 8.23306s + .9590716 = 0. 

When z has been computed, x is found from 

2000s + 80 
100s — 1 


x 
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Only linear transformations which make the coefficient of 
the second term of the complete equation or of the equation 
reversed zero or small are used, as other transformations 
become too complicated to make the method practicable. 

IV.— Transcendental Trinomial Equations. 

Let an equation of the form y + af(y)+b = 0, where / (y) is 
a transcendental function, be called a transcendental trinom¬ 
ial equation. Such equations are readily solved by the 
method, provided the resulting series is rapidly convergent, 
but in the absence of a transformation which insures rapid 
convergence the method has little practical value. 

Suppose the equation 2y + \ogy — 1000 = 0 to be given. 
Applying the method to 2?/ + .Tlog?/ — 1000 = 0, if the 

Napierian logarithm of y is taken, ?/ 0 = 5000, d ~ = — 4.30625, 

tlx* 


i 7^2 = +0.000215, and ?/ = 4995.G9; if the common logarithm 

ttXft 

of y is taken, 2 / 0 = 5000, + = —l.S494S,*^ = + 0.00018, and 

UX 0 - u.T 0 “ 

?/ = 499S.15. 


V.— Expansions. 


If y=z+v<f(y), where v and z are independent variables, 
Lagrange’s series expands any function of y in powers of v. 
These expansions may be obtained by ;writing y=z+vx<p(y ) 
and expanding / (y) y which now becomes a function of x, by 
Maclaurin’s series and making x unity in the result. 

The method will be illustrated by obtaining two expan¬ 
sions which occur in theoretical astronomy. From the 
equation E = M + e sin E , where E is the eccentric anomaly, 
M the mean anomaly and e the eccentricity of the orbit, it is 
necessary to find E and (1 — e cos E) 2 . 

To find E y write E = M + cx sin E , whereby E becomes an 
implicit function of x. Differentiating twice in succession 
with respect to x, 


dE • r, , jpdE 

~r~ = e sin E + ex cos E 1 , 
dx dx 

d 2 E dE d-E 

= 2 e cos E _— f- ex cos E -■ — ex 

dx 2 dx dx 1 


ex s 
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Making x zero, E 0 =M, , 0 = e sin M, = 2c 2 cos .1/ sin M. 

UXq (IXq 

Substituting in Maclaurin’s series and making x unity, 

E — 3f + e sin 3f -f sin (2Af) -f- . . . . 

To find (1 — ecos E)~ 2 , write E = M+ex sin E and 
y = (1 — e cos £)~ 2 . Since y is a function of x through E y 

— 2e (1 — e cos E )- 3 sin 

= Ge 2 (1 — e cos sin 2 E c -~ 
ax 1 eta 

— 2« (1 —geos#)—3cos #(^|) 2 

— 2c (1 — e cos ^7)— 3 sin E ~j- 

Placing x = (), when E = M, = e sin M and 

< ^=2e 2 sin M cos M, 
dx 2 

Vo — (1 — 6 cos J/)- 2 , 


= — 2e 2 (1 — e cos J/)— 3 sin 2 J/, 

= G^ 4 (1 — e cos 31 )— 4 sin 4 31, 

dx o 2 

— 6c 3 (1 — e cos J/)— 3 sin 2 31 cos 31. 

Substituting in Maclaurin’s series and making x unity, 

(1 — e cos E)~ 2 = (1 — e cos 31)— 2 — 2e 2 (1 — e cos J/)— 3 sin 2 J/ 
+ 3e* (t — e cos 31 )—* sin 4 31 
— (1 — e cos J/)— 3 sin 2 J/ cos J/ + . . . . 

In like manner all expansions obtained by Lagrange’s 
series maybe obtained by a direct application of Maclaurin’s 
series. Of course it is evident that if e is considered a varia¬ 
ble the derivatives with respect to e may be formed and the 
introduction of x is unnecessary. 
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Historical Note. 

Lagrange, in the memoir “ Nouvelle methode pour 
resoudre les Equations Litterales par le moyen des Series,” 
read before the Berlin Academy in 1770, found all the roots 
of an equation in infinite series. McClintock, in Volume 
xvii of the American Journal of Mathe?natics , obtained by his 
Calculus of Enlargement series better adapted to computa¬ 
tion. It was recognized that these series may be obtained 
by Lagrange’s series. McClintock calls the coefficients of 
the terms which have been underscored the dominants of the 
equation. The method of the present paper brings the com¬ 
putation of the roots of equations by means of series within 
the range of elementary instruction. 

Since completing this paper the author found in an extract 
of a letter from Cauchy to Coriolis, of January 29, 1837, 
published in the Comptes Rendus of the Paris Academy, an 
announcement of important results to be obtained by break¬ 
ing up an equation into two parts and introducing as a factor 
a parameter into one part, which parameter is ultimately to 
be made unity. In a postscript Cauchy states he discovered 
the advantage of making one part a binomial. But the 
author has been unable to find the method sketched in this 
letter developed. It would indeed be surprising if a method 
so strikingly direct had escaped notice. 

Lehigh University, South Bethlehem, Pa. 


